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Two-level quantum systems, qubits, are uot the only basis for quantum computation. Advantages 
exist in using qudits, d-level quantum systems, as the basic carrier of quantum information. We show 
that color codes—a class of topological quantum codes with remarkable transversality properties—can 
be generalized to the qudit paradigm. In recent developments it was found that in three spatial 
dimensions a qubit color code cau support a transversal non-Clifford gate, and that in higher spatial 
dimensions additional non-Clifford gates can be found, saturating Bravyi and Konig’s bound [Phys. 

Rev. Lett. 110 , 170503 (2013)]. Furthermore, by using gauge fixing techniques, an effective set of 
Clifford gates can be achieved, removing the need for state distillation. We show that the qudit color 
code can support the qudit analogues of these gates, and show that in higher spatial dimensions a 
color code can support a phase gate from higher levels of the Clifford hierarchy which can be proven 
to saturate Bravyi and Konig’s bound in all but a finite number of special cases. The methodology 
used is a generalisation of Bravyi and Haah’s method of triorthogonal matrices [Phys. Rev. A 86 
052329 (2012)], which may be of independent interest. For completeness, we show explicitly that the 
qudit color codes generalize to gauge color codes, and share the many of the favorable properties of 
their qubit counterparts. 


I. INTRODUCTION 

Quantum technologies are often developed in the qubit 
paradigm, where the basic carrier of quantum information 
is a two-level quantum system. Qubits are a natural choice 
because binary is the language of classical technologies. 
However, even here, despite the prevalence of binary, its 
supremacy is questionable. Indeed, Donald Knuth has 
advocated the use of balanced ternary, a 3-state classical 
logic [1] . In the quantum domain, qudits offer a state space 
with a richer structure than their two-level counterparts, 
and the merits of this for quantum information have been 
explored in many contexts HHin]. Recent experiments 
have shown even large d quantum systems can be precisely 
controlled mm- 

Here we present a qudit generalisation of a power¬ 
ful class of quantum error correcting codes, the color 
codes [ISUTT] . Along with surface codes [HI [19] , they 
constitute the most successful topological codes. Us¬ 
ing topology, quantum codes have achieved high error 
thresholds, whilst proving more practical than concate¬ 
nated codes. Thresholds of qudit surface codes indicate 
improvements with qudit dimension (TH). Qubit color 
codes have several advantages over qubit surface codes, 
and we show these features can be transferred over into 
the qudit setting. 

There has been some prior investigation into qudit color 
codes for prime-power d (SU] and color codes based on 
more general groups |21] . although these works were re¬ 
stricted to 2D topologies. In this paper, we generalize 
color codes to any qudit dimension d and spatial dimen- 
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sions /i up to the point where pt factorial is a multiple of 
d. Specifically, given any lattice suitable for constructing 
qubit color codes, we show how to use the same lattice to 
construct a qudit color code. For qubits, a non-Clifford 
gate can be implemented in color codes in 3 and higher 
spatial dimensions transversally, i.e. by a tensor product 
of local unitary gates, an inherently fault-tolerant pro¬ 
cedure To avoid confusion between the spatial 

dimension of the lattice and the Hilbert space dimension 
of the qudit, we shall always denote the former by the 
letter pi and the latter by the letter d. 

Recently, it has been shown by Bravyi and Konig [5S] 
that a quantum error correcting code in pt spatial dimen¬ 
sions can support a gate with constant depth from at 
most the /rth level of the Clifford hierarchy. The fact 
that color codes can be shown to saturate this bound 
with transversal gates is a very promising feature, and 
when combined with gauge fixing techniques [miiMj 
enables universal quantum computation without the need 
of magic state distillation (5] |6l [29l[34| . The structure of 
the qudit Clifford group is very different from its qubit 
counterpart 11135 ]. Nevertheless, we find that 3D color 
codes also provide transversal non-Clifford gates in the 
qudit case. 

Recently the color codes were generalized to gauge color- 
codes m —subsystem codes with many advantageous 
features—including low weight error detection measure¬ 
ments, universal transversal gates via gauge fixing for 
pL > 2, fault-tolerant conversion m between codes of 
different spatial dimension [ 35 ], and for the pL = 3 case, 
single shot error correction m — a robustness to measure¬ 
ment errors without the need for repeated measurements. 
We show that the qudit color codes introduced here can 
also be generalized to gauge color codes. 

The main technique which we employ is a bipartition of 
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the vertices in the graph that defines the code into starred 
and unstarred vertices. We call this the star-hipartition, 
to distinguish it from the other important colorings which 
define the color codes. The commutation properties of 
the stabilizer and logical operators of the color codes (and 
gauge color codes) in the qubit setting can be reduced 
to the fact that the pairs of operators X X and Z ® Z 
commute. The star-bipartition we introduce replaces 
some operators with their complex conjugate with respect 
to the computational basis of Z eigenstates. For example, 
replacing the above operators with X 0 X* and Z Z) Z*, 
respectively. Crucially, this latter pair of operators does 
commute for qudits of any dimension, and this becomes 
the starting point for a generalisation of the color codes 
from qubits to qudits of any Hilbert space dimension d 
and any spatial dimension p. The sharp-eyed reader will 
note that X* = X. We write the star explicitly on a real 
matrix here, and throughout, to emphasize symmetry and 
to simplify notation. Thus we see that as long as a pair 
of X- and Z-type stabilizers have in common an equal 
number of starred qudits, they will commute. 


Furthermore, the star-bipartition provides a general 
framework for constructing transversal gates from higher 
levels of the Clifford hierarchy. While elements of this 
technique can be seen in earlier work nnno], this is 
the first time that it has been exploited systematically. 
The second key technical component of our work is a 
generalisation of the triorthogonal matrix technique by 
Bravyi and Haah m- 


This paper is structured as follows. We start in Sec. [IT] 
by reviewing the stabilizer formalism for d-level systems. 
In Sec. Ill we describe how to generalize a qubit color 
code in arbitrary spatial dimensions to a qudit color 
code by employing the notion of star-conjugation. In 
Sec.jTVjwe review and generalize the triorthogonal matrix 
construction. In Sec.jVjwe derive conditions on the lattice 
that must hold for a transversal non-Clifford gate to 
be implementable on the code. Sec. |VI| explains how 
the Hadamard can be implemented on the same lattice, 
although not in the color code, using the technique of 
gauge fixing. In Sec. |VII| we show gauge color codes can 
naturally be defined for all the codes we have introduced, 
inheriting the favorable features of the qubit gauge color 
codes. We conclude in Sec. Eml 


II. QUDIT STABILIZER CODES AND THE 
CLIFFORD HIERARCHY 


The qudit analogues are 

^ = XI1-^' + (ji: ^ = X 1-^'^ ’ 

ieZd iGZd 

( 1 ) 

where w = e^. With a slight abuse of terminology, 
we shall say that two operators A and B “w-commute” 
if AB = uiBA and note that w-commutation holds for 
X and Z. These generalized operators simplify to the 
familiar Pauli operators for d = 2. 

The single qudit Pauli group is generated (up to global 
phases) by X and Z and the n-qudit Pauli group P®" 
is the n-fold tensor product of the single qudit Pauli 
group. Consider an abelian subgroup S C 7^®", such that 
uj^l ^ S for nonzero j, then we say that 5 is a stabilizer 
group, and we refer to its elements as the stabilizers. 
The stabilizer group defines an error correcting code with 
codewords corresponding to states \ip) that are stabilized 
by the stabilizers, i.e. S \ijj) = \tjj) for all S G S. 

The logical operators correspond to the set of operators 
that commute with S but are not contained in it. A pair 
of logical operators Xi and Zi w-commute with each other 
and hence encode one qudit. 

Gottesman and Chuang [35] introduced a classification 
of quantum gates known as the Clifford hierarchy (CH), 
which can be defined recursively [39] as 

Vi = {C/|P^t7PC/1' G Vi-i V P G Pi}, 

( 2 ) 

where Pi is the Pauli group. For example, P 2 is the group 
of operators that leave the Pauli group invariant under 
conjugation and is called the Clifford group. We shall 
describe some important Clifford group gates H, S and 
A(A) below. In prime qudit dimensions, these are known 
to generate the whole Clifford group [H HO] ■ The gate H 
is the qudit version of the Hadamard gate (also known as 
the discrete Fourier transform), 

X (3) 

''“i.fceZd 

the S gate is the generalisation of the qubit 7r/4-phase 
gate, 

^ = X|j) 01. (4) 

je^d 

and A(A) is the controlled-A gate (also known as the 
SUM gate), 


We will consider d-level quantum systems (qudits) as 
the building blocks for the constructions of quantum 
codes |3]. Unless stated otherwise, the qudit dimension 
d is assumed to be any integer greater than two. The 
conventional basis states for the d-level system are taken 
to be |j), for j G The single qubit Pauli matrices X 
and Z have natural extensions in higher dimensions |3]. 


A(^)= X \j)c\k(Bl)Aj\c{k\t^ (5) 

j,keZd 

where c and t are the control and target qudits, respec¬ 
tively. 

The set of gates Vi in the hierarchy contains all gates 
from lower levels of the hierarchy. To refer to gates in 







3 


level I of the hierarchy but not level / — 1 we shall say 
that the level I of the hierarchy is the lowest level of the 
hierarchy for which the gate is a member. 

In prime dimensions, it is known that it suffices to 
supplement the Clifford group with just one non-Clifford 
gate from the third level of the CH in order to obtain 
a universal set of gates [5]. Such a gate is not unique, 
and in the qubit case, the T gate diag(l, is usually 

chosen for this purpose. For the qudit case, we choose 
the following particularly convenient definition for the 
T gate, which is valid in all dimensions except when 
d = 2,3,6 HIH], 

T = ^ \j) {j\ . (6) 

It is a consequence of lemmas and (below) that for 
d ^ 2, 3, 6 this gate is non-Clifford and inhabits the third 
level of the CH. In d = 3 and d = 6, the gate defined 
in equation ^ is not non-Clifford, as it reduces to the 
Pauli Z gate since = j mod 3 and mod 6. However, 
the following definition provides a suitable alternative T 
gate for these dimensions. The gate is non-Clifford and 
in the third level of the CH: 

73,6 = 7 ^' 1 -^'^ 

where = uj and where the function in the exponent 
is evaluated in regular arithmetic (or equivalently modulo 
3d). When we refer to “the T gate” in this paper we will 
always mean a gate of the form of equations ^ or 0, 
depending on the qudit dimension under consideration. 
For notational convenience we suppress the dependence 
of T on d, since it will always be clear by the context 
which T gate we require. 

Notice how the T gate has a cubic power in the exponent 
of w, in contrast to the quadratic power in the case of 
the S gate. For prime dimensions, the first investigation 
characterising all the phase gates from the third level 
of the CH was performed by Howard and Vala [3]. In 
general, there is a close correspondence between the order 
of the polynomial in the exponent of w and the lowest level 
of the CH the phase gate belongs to. Let us define the 
following family of phase gates in terms of a polynomial 
function fr{j) of degree r, such that r < d with coefficients 
am G ^d, so frU) ■= TJm=0<^rnj”"- 

^.;/^«|j)(j|, (8) 

je^d 

one can then prove the following useful lemmas. 

Lemma 1. For all d, all r < d and all functions fr{j), 
the gate Rf^ is in the rth level of the Clifford hierarchy. 

Proof. The proof of this is simple and concise. We begin 
by calculating 

Rf^_, = (9) 


where/,._i is a new function/r_i(j) = /r(j + l)-/r(j) = 
rarj'~~^ + • ■ ■ ■ Expanding out, the degree r terms cancel, 
so the leading term is rarj^~^ and extra terms are all 
degree r — 2 or smaller. We now observe that if r = 
1, Rf^ is a Pauli operator. Using the definition of the 
Clifford hierarchy in equation 0, the lemma follows by 
induction. □ 

Lemma 2. For all d, all r < d and all functions fr{j) 
satisfying rlor ^ 0 (mod d), the gate Rf^ is not in the 
(r — l)th level of the Clifford hierarchy. 

Proof. If the gate is in the (r — l)th level of the hierarchy, 
then applying the inductive transformation in the previous 
proof r — 1 times must return a Pauli operator, and thus 
applying it r times results in an operator proportional 
to the identity. By chaining together transformations 
of the form of equation 0 n times, we obtain /i(j) = 
rlorj + c. If the corresponding operator is proportional 
to the identity then this function must be constant, and 
thus rlor = 0. □ 

These two lemmas provide us with a simple way to 
generate gates at all levels of the Clifford hierarchy as 
required. They fail when rlor = 0, which was the case 
above for r = 3 and d = 2, 3,6. In those exceptional cases, 
gates can be discovered by moving to higher roots of unity 
as illustrated by equation 0. 

III. QUDIT COLOR CODES 

Color codes [13] are a class of topological qubit stabilizer 
codes that be defined on a topological space of any spatial 
dimension |42j /r > 2. The /r-dimensional manifold is 
celluated into a lattice of objects called k-cells for all 
spatial dimensions 0 < k < ^. For example, vertices are 
0-cells, edges are 1-cells and a 2-cell is a cell defined on 
the faces, or plaquettes, of the lattice. We define 

Definition 1. A lattice C is called a yi-colex (colex is 
short for “color complex”) whenever 

1. it is a celluation of an orientable y,-dimensional 
manifold without a boundary; and 

2. every vertex has {y + 1) neighbors ({y + 1)-valency); 
and 

3. the y-cells are {y -\- 1)-colorable. 

Any /r-colex defines a qubit color code. For example, 
the smallest 3-dimensional color code is a 15-qubit code 
defined on the lattice illustrated in figure 0a). The 
stabilizer group is generated by face operators (with a 
Z operators assigned to the vertices around the face or 
2-cell) and cell operators (with an X operators assigned 
to the vertices around each 3-cell), see figure 0b). Our 
results show a ^-colex also defines a qudit color code. An 
alternative description of a /i-colex is that its dual C* is 
a simplical lattice. This is lattice where each fc-cell is a 
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simplex, an object with fc + 1 vertices. In the dual, we 
have conditions: (2*) every ^-cell has + 1) neighbors; 
(3*) the vertices are {fi + 1) colorable. 

Such codes can be constructed m by starting with a 
closed hyperspherical lattice and then removing a vertex 
to “puncture” the surface. Alternatively, as we show in 
Sec. lmcl one can construct them directly by defining a 
suitable boundary on a regular lattice structure. Such a 
code encodes a single qubit and has the attractive feature 
of a transversal non-Clifford T gate. The qudit color 
codes have many useful properties. For more details, we 
refer the reader to 

The stabilizer generators are guaranteed to commute 
in this construction. Those of the same type {X or Z) 
trivially commute, while those of different types commute 
since they always meet in pairs. As already remarked, 
X Z) X and Z Z) Z commute only in d = 2. However, 
by replacing one operator of each pair with its complex 
conjugate we produce a pair of operators that commute 
for all d. These are X (g) X* and Z Z*. To define color 
codes in all qudit dimensions, we need, therefore, to find 
a construction which allows us to take advantage of the 
commutation of A(g)A* and Z®Z*. This can be achieved 
by identifying and exploiting a bipartition or bicoloring 
of the graph defining the lattice. To avoid confusion with 
other colorings of the lattice important for color codes we 
call this the star-bipartition. 

We note that a similar construction was defined (for 
prime power qudit dimension only) in [20j . Here we go 
further and show that the star-bipartition is the starting 
point for an identification of a broad family of transversal 
gates on these codes, including non-Clifford gates satu¬ 
rating Bravyi and Konig’s bound. 

A. Star-bipartition 

The star-bipartition is a bipartition of the color code 
lattice C. It thus divides the set of lattice vertices into 
starred and unstarred vertices where neighboring vertices 
always belong to different sets. Here we prove the follow¬ 
ing useful bipartition lemmas 

Lemma 3 (Star-bipartition lemma). Let C he a fj,-colex, 
then its vertices can he 2-colored into starred and unstarred 
sets, u* and v,, respectively. 

Recall that the /x-colexes have several properties listed 
in Def. that are essential to proving the above lemma. 
Using I ... I to denote the number of elements in a set, we 
also have 

Lemma 4 (Starring of cells lemma). Let C he a pL-colex, 
so that its vertices are partitioned into u, and u* according 
to lemma^ It follows that 

1. |u*| = |u,|; and 

2. any k-cellC with 0 < k < p contains a equal number 
of vertices from each partition, so jCnu*! = |Cnu, |. 


Lemmas and or variants thereof, were already 
proved in prior research na [ig na |42]. These results 
are especially important to qudit color codes and play a 
fundamental role, and so, for completeness, we present 
our own proofs. Our presentation of lemma ^ has the 
merit of being more concise than that of Ref. Q. 

Before presenting our proof, we review the concept of 
orientation from topology theory [43) . An orientation on 
a /x-simplex is an ordered list of its vertices, with orien¬ 
tations considered equivalent whenever they differ by a 
permutation generated from an even number of trans¬ 
positions (a transposition is swap of two elements). All 
permutations can be generated by either an even or odd 
number of transpositions, and so there are two possible 
orientations for a ^-simplex. The /x-simplex contains sub- 
simplices that can be obtained by removing a single vertex. 
This vertex removal provides an induced orientation on the 
subsimplex, such that removing Vj from a simplex with ori¬ 
entation {vi,V 2 ,..., Vj-i,Vj,Vj+i ,..., Vfj^, u^+i} induces 
an orientation {vi,V2, ■ ■ ■, Uj_i, ... ,v^ , u^+i}. Con¬ 
sider two oriented /x-simplices with a common subsimplex 
(that is, they are are neighbors). We say their orientations 
are consistent if they induce opposite orientations to their 
common subsimplex. An entire lattice of simplices is said 
to be orientable if there exists a choice of orientations 
such that all pairs of simplices are consistently orientated. 
An example of a consistently orientated lattice is shown 
in Fig. Orientablity of a lattice is a topological feature 
that depends on the underlying manifold and from this 
definition one can show that many familiar manifolds are 
orientable, including Euclidean, spherical and toroidal 
manifolds. We now employ these concepts in our proof. 

Proof. To prove the star-bipartition lemma, we switch to 
the dual lattice C*. In the dual picture, we need to show 
that the /x-dimensional cells can be 2-colored. Property 2 
of Def. [T] ensures that the dual lattice is simplical. Prop¬ 
erty 3 of Def. [^ensures the dual vertices are (/x-l-l) colored. 
This coloring can be described by a map C : £* —>■ 
that assigns color label C{vj) to vertex Vj. Given a sim¬ 
plex (T with an orientation, say {vi, V 2 ,. ■. we can 

apply the color map to the orientation to obtain a colored- 
orientation {C{vi),C{v2), .. .C(x;^+i)}. This defines a bi¬ 
partition (the star-bipartition) with two simplicies belong¬ 
ing to the same partition if their colored-orientation differs 
by an even permutation. The remainder of the proof es¬ 
tablishes that neighbouring simplicies belong to different 
partitions. Two neighbouring oriented simplices A and B 
always share a common (/x — 1) subsimplex. Let us call 
the vertices they share in common {vi, ..., v^}. Simplex 
A also contains a vertex va and we label its orientation 
{vi, ... ,Vfj,,VA}. Similarly, simplex B contains additional 
vertex vb and has orientation {un(i), x'n( 2 )) • ■ • ^n(Ai )5 
where H is some permutation of vertex labels. No gen¬ 
erality is lost by placing va and vb last in their respec¬ 
tive orientations as this can always be achieved with an 
even number of transpositions (assuming p > 1). Since 
the manifold is orientable, we can assume a consistent 
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FIG. 1: (Color online) A patch of an orientated simplicial 
lattice with 3 colored vertices. The orientation is visually 
represented by a clockwise symbol. Formally, the orientation 
is an ordered list of vertices {x, y, z] moving clockwise round 
the simplex. For neighboring simplices A and B, we also show 
the induced orientation on their edges. For an edge orientation 
{x, y}, we illustrate this with an arrow going from x to y. The 
simplices A and B induce opposite orientations to their com¬ 
mon edge, and so are consistently orientated. If a simplex has 
an orientation {a, b, c} with colors {red, green, blue} then we 
label it starred, and otherwise unstarred. The star-bipartition 
lemma further formalises this notion of starring in terms of 
color and orientation. 


orientation on their common subsimplex. We see that 
A induces the orientation (rii,..., and B induces 
the orientation {t^n(i)) • • ■ j t'n(;i)}j and so consistency de¬ 
mands that n is an odd permutation of the vertex la¬ 
bels. Next, we show this lifts to odd permutation in 
the color-orientation of A and B. Simplex A has color- 
orientation {C(fi),... ,C{Vfj,),C{vA)} and simplex B has 
color-orient at ion {C(t>n(i)), ■ • • Since 

the lattice is (/i + I) colored, we know va and vb must 
possess the same color, which is whatever color is absent 
from their common subsimplex. Permutation of vertex 
labels results in a corresponding permutation of colors. 
Therefore, the color-orientations of A and B differ by an 
odd permutation of their first fj, elements, which corre¬ 
sponds to an overall odd permutation. Since all neigh¬ 
boring /i-simplicies must have opposite color-orientations, 
they belong to separate partitions of our star bipartition, 
proving the lemma. □ 

Let us now turn to lemma |4l 


Consider any fc-cell C. It defines a sublattice Cc with 
same vertex set as C. If Cc is a fc-valent with k > 0, 
then the above argument applies and C contains equal 
number of starred and unstarred vertices. The following 
reasoning parallels that in Ref. m- We show explicitly 
that Cc is a fc-colex for the k = yc — 1 case, but any k 
is reached by iteratively applying the argument down to 
the desired k. Again, we must switch to the dual lattice 
where C*q is obtained by removing a single vertex C* from 
C*. The sublattice retains all simplices containing C *, 
but with their dimension reduced by the removal of C*. 
Therefore Cq is a simplical lattice of dimension (p — 1) 
and so C is /i-valent. Our proof only requires the correct 
valency, which we have shown, but note that Cc also has 
the correct coloring for a (/r — l)-colex. Specifically, if 
C and C intersect on some (/r — l)-cell of Cc, then we 
assign it the color of C". □ 

The above results concern a lattice without a boundary. 
However, if one punctures the code removing a starred 
vertex then we have 

Corollary 1. Let C be puncturing of a yi-colex, with the 
inherited bipartition into u* and v,. Then 

1. |z;*| = |t^,| - 1; and 

2. any k-cell C with 0 < k < y, contains a equal number 
of vertices from each partition, so |C'nz;*| = ICnt),!. 

Property 1 is immediately follows from point 1 of 
lemma as a single starred vertex has been removed. 
Property 2 is identical to its partner in lemma The 
punctured lattice only keeps cells that did not contain the 
punctured vertex, and so the property is directly inherited. 
The dual lattice C* has been an essential proof tool for 
these lemmas, but in the rest of the paper we shall only 
consider the primal lattice. 


B. Qudit color codes in arbitrary spatial dimension 

Before we write down the stabilizer generators for these 
codes, we shall use the star-bipartition to define notation 
for an important family of transversal operators which 
we call star-conjugate transversal. Throughout this pa¬ 
per, when discussing color codes we identify qudits with 
vertices in the lattice defining the code. Let v denote a 
set of vertices (and thus the corresponding qudits). Let 
us introduce the notation 


Proof. We shall begin with part 1 of the lemma. Consider 
a bipartite lattice that is r-valent. We count the number 
of edges Ne- Every edge is incident to one, and only one, 
starred vertex. Furthermore, each vertex is contained 
in r edges, and so the total count is Ne = |r’*|r. The 
same argument applies to unstarred vertex, so Ne = |u, |r. 
Equating |ri*|r = It'.lr, we see that r > 0 entails |z;*| = 
jf;, |. Since a y-co\ex is {y -I- l)-valent, we have proven 
part 1 of lemma 


U[v] = (10) 

to denote the tensor product of unitary U acting on each 
qudit identified with the vertices in v. 

We call an operator star-conjugate transversal when it 
has the following form 

U = U[v,]®U*[v4. 


( 11 ) 
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FIG. 2: (Color online) (a) The smallest instance of the qudit 
3D color code with the X-type stabilizers colored red, green, 
blue and yellow. The 1-cells (edges) of the code are also colored. 
The vertices are also colored so the set v, is represented by 
black circles and the set n* is represented by white stars, (b) 
A single Z stabilizer of the tetrahedral 3D color code. The 
plaquette can take the color yellow if considered as a face of 
the green 3-cell, and vice versa, (c) A single X stabilizer of 
the tetrahedral 3D color code. 


and unstarred vertices, which is a point discussed further 


in Sec. IV C Hence these stabilizer generators commute 
as required. Note that the qubit d = 2 case is included in 
our definition. 

We noted above that, in the constructions we consider, 
the code will contain one more unstarred qudit than there 
are starred qudits. We can thus define the logical encoded 
Pauli operators for the code star-conjugate transversally. 


V = V = V[n.](g)V*K] 

Z = Z = Z[vt\ ® Z*[vi,]. 


(14) 


The fact that |u*| = |z;,| — 1 ensures that the logical 
operators satisfy the same commutation properties as X 
and Z. One can verify that these operators commute with 
the stabilizer operators defined above by recalling that 
each 2-cell P and and /i-cell C contains an equal number 
of starred and unstarred vertices. 


C. Constructing codes of any code distance 


In other words, U consists of U applied to all unstarred 
vertices and U* applied to all starred vertices. 

To define qudit color codes we need to introduce two 
types of stabilizer generators which are defined with re¬ 
spect to cells of different dimensions within the lattice. 
As in the qubit case, we associate .^-type stabilizer gener¬ 
ators with the ^'-cells of the lattice and V-type stabilizers 
with the {fji — fj! + 2)-cells [22]. In this section, for sim¬ 
plicity of presentation, we will take the example of /i' = 2 
so the Z-type stabilizers are associated with plaquettes 
(2-cells) and V-type stabilizers are associated with /i-cells, 
although the construction holds in the more general case, 
which we shall consider later. For example in 3D, the X 
stabilizers act on the vertices contained in a 3-cell of the 
lattice, and the Z stabilizers act on the vertices contained 
in a 2-cell (plaquette). This can be seen in figure]^ 

Setting fi' = 2 the stabilizer generators therefore take 
the form 


Sx,c = x[u.nc]®v*KnC'], (12) 

Sz,p = z[v,np](z z*[v^np]. (13) 

for all 2-cells P, and /i-cells C of the lattice. Recall that 
we write the conjugate operator X* explicitly to emphasis 
the ubiquity of starring, even though AT is a real operator 
in the computational basis and X = X*. It is not just the 
Sz,p operators where starring is non-trivial, but many 
logical operators also require starring. 

Let us denote the group generated by Sx,c operators 
Sx and Sz,p operators by Sz- All elements of Sx com¬ 
mute as do all elements of Sz- It remains to show that all 
elements of Sx commute with all elements of Sz- This 
follows from the fact that X 0 X* commutes with Z Z) Z*. 
The above construction ensures that whenever cell C and 
cell P overlap they overlap on an equal number of starred 


In topological stabilizer codes we expect to be able to 
increase the number of physical systems encoding the 
quantum information in order to protect the information 
more effectively—a property characterized by the code 
distance. We have remarked already that suitable codes 
can be constructed by puncturing a hypersphere. In this 
section we provide, however, a constructive alternative 
method to show that in /i spatial dimensions the color 
code distance can be made arbitrarily large. 

The code distance is the weight of the smallest logical 
operator. It was convenient above, in equations (14), to 
define the logical operator basis star-conjugate transver¬ 
sally. However, multiplying by a subset of the stabilizer 
group we can localize Z to any edge of the polytope formed 
by the lattice, where such an object is a 1-dimensional 
sub-manifold of the lattice between two of the vertices 
of the polytope. Each such vertex is contained in only 
a single /r-cell (as opposed to those contained in two or 
more cells of the lattice), and hence there are exactly 
(/r -I- 1) vertices of the polytope. 

In /r spatial dimensions, an edge of the polytope is 
comprised of 1-cells of two different colors. Similarly, 
2-cells contain 1-cells of two colors. Beginning from the 
transversal definition of the operators in equations (14) 
and multiplying Z by the subset of Z stabilizers defined 
on 2-cells containing 1-cells of the same pair of colors, the 
logical operator is localised to the edge of the polytope 
uniquely defined as containing 1-cells matching this pair of 
colors. The constraints placed on the subset of stabilizers 
chosen, and the initial choice of transversal operator basis 
ensures that the minimum length of a string-like logical 
operator defined in this way is the length of the edge of 
the polytope colored as indicated. The minimum length 
of a string-like Z in such a definition, and hence the 
code distance, is then length of the shortest edge of the 
polytope. Thus by showing that the length of an edge of 
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FIG. 3: (Color online) To increase the distance of the code, we 
tile 3D space with truncated octahedron cells (with appropriate 
boundary conditions). Here we show a portion of the lattice. 
A red cell fits into the hollow formed by the green, blue and 
yellow cells towards the bottom of the figure, and thus the 
tiling proceeds, ensuring the 4-colorability of the lattice is 
preserved. 


the polytope can be made arbitrarily large we demonstrate 
the ability to create codes of any distance. 

We begin with an example in 3D and then generalize 
the argument to higher spatial dimensions. The particular 
choice we make for this example is a 4-colorable tiling of 
truncated octahedra, part of which is illustrated in figure!^ 
This is not a unique choice—for an alternative see |4^ 
where the 3D color code is defined on a lattice constructed 
from cubes and truncated octahedra. 

In addition to the two lattices already mentioned, we 
have identified two further possible choices of regular 
tilings in 3D Euclidean space on which a color code may be 
defined. One is the cantitruncated cubic honeycomb, com¬ 
prising truncated cuboctahedra (polyhedra with square, 
hexagonal and octahedral faces), truncated octahedra, 
and cubes. The second is the omnitruncated cubic honey¬ 
comb comprising truncated cuboctahedra and octahedral 
prisms. This list is not exhaustive, for instance many 
more tilings may exist in curved space. 

The method proceeds by carefully cutting a block of 
the lattice out, in order to form a tetrahedron, similar 
to that illustrated in figure [^a). The main difference is 
that the new tetrahedron contains more qudits and can 
have an arbitrarily large (though always odd, in order 
to ensure the w-commutation of the logical operators) 
number of qudits contained in its edges. The shape cut 
from the lattice may not look tetrahedral at first but 
may be deformed to a regular tetrahedron, with vertices 
belonging only to a single cell forming the four vertices 
of the tetrahedron. 

In ^ spatial dimensions a suitable tiling of the space 
as outlined at the beginning of section |III| must be found. 
Besides these requirements placed on the lattice construc¬ 


tion, the “block” (polytope) cut out of the lattice to form 
the higher-dimensional analogue of the tetrahedron must 
adhere to certain rules. 

1. The polytope formed must have (/r-|- 1) boundaries, 
each of a different color. The color of the boundary 
corresponds to the color of the string-like operators 
that can end on that boundary. 

2. Every edge of the polytope must contain an odd 
number of qudits. As stated above, this is to en¬ 
sure the w-commutation of the logical operators is 
preserved. 

3. There should be exactly /r vertices that belong only 
to a single cell of the lattice. These cells should 
be different colors and these vertices will form the 
vertices of the polytope. 


D. Error detection 

Errors are detected by measuring all the stabilizer gen¬ 
erators. We shall not present a detailed analysis here, 
since errors and syndrome are related for these codes 
in a similar way to the qubit case. For example, in 3 
spatial dimensions a single Z-type error on an unstarred 
qudit will result in a measured eigenvalue of w in the X 
stabilizers corresponding to the four cells containing that 
vertex (see figure]^. Similarly, a single A-type error on 
an unstarred qudit will result in a measured eigenvalue 
of w in the six faces that contain the vertex (see figure]^. 
In general, a or error will lead to a measurement 
outcome of Syndromes arising from sets of multiple 
errors can be calculated using standard techniques. 

A classical algorithm called a decoder may be used to 
interpret the syndrome and infer the correction operator 
that must be applied to return the code to its original state. 
There are many examples of decoding algorithms for 2D 
qubit color codes that have been developed recently m- 
I50j , but very little is known for the case of higher qudit 
and spatial dimensions. 

However, in the case of surface codes, it has been shown 
that some proposed renormalization group (RG) algo¬ 
rithms can be generalized to qudit codes in a straightfor¬ 
ward manner IS1IS2]- Moreover, recent work for the 3D 
fault-tolerant implementation of the qudit surface code 
using a hard-decision RG decoder [5] suggests that adap¬ 
tations of such algorithms for higher spatial dimensions 
should be possible, however, the error correction thresh¬ 
olds would likely be degraded as the spatial dimension 
increases due to the larger stabilizer generators. 

We shall not present a decoder in this paper, but re¬ 
mark that, unlike in qubit codes where neighboring errors 
lead to cancellation of syndromes in between, such cancel¬ 
lation becomes unlikely with increasing qudit dimension, 
and thus in general more information is available to the 
decoder. In fact, it has been shown that it is possible to 



design decoders to exploit this additional information to 
obtain higher thresholds with respect to qubit codes [9]. 

For the case of gauge color codes, very recently a hard- 
decision RG decoder m was implemented for the qubit 
3D construction. As with the color code, we expect a 
generalization of such an algorithm to gauge color codes 
in arbitrary spatial and qudit dimensions to be straight¬ 
forward. 


IV. m*-ORTHOGONAL MATRICES AND 
CODES 

Brayvi and Haah introduced a powerful framework 
for defining codes with a transversal non-Clifford gate; 
codes defined in terms of triorthogonal matrices. Here we 
introduce a significant generalisation to Bravyi and Haah’s 
approach, extending to qudits codes with transversal gates 
from higher up the Clifford hierarchy. 

A. Matrix representation of quantum codes 

We begin by reviewing how Bravyi and Haah represent 
CSS codes with matrices m- They use a matrix G over 
Z 2 that has linearly independent rows under modulo 2 
arithmetic. The matrix G is broken up into two blocks 
Gi (with k rows) and Gq (with s rows). This dehnes 
a quantum code with k logical qubits. The elementary 
logical basis state | 0 l) is written 

|0l) = X! 

/Gspan[Go] 

= i E (15) 

where denotes matrix multiplication modulo d, and we 
take advantage of the fact that the rows of Gq are linearly 
independent to represent the terms in the superposition 
by a row vector y transposed and multiplied with Gq. For 
the other logical computational basis states la;^) where 
a: G Z 2 , we have 

\xl) = ^ ^ ly'^-GoOx'^-Gi), 

(16) 

where addition of row vectors is elementwise modulo 2 as 
acknowledged by the symbol 0. Note that in the special 
case that k = 1, x is a. scalar (or a 1 x 1 row vector) and 
the transpose operation is trivial. For qudits, quantum 
codes are defined by a matrix G taking elements from Z^, 
so that 

\ xl ) = 0= V IG ■ Go ® ■ Gi). 


The key differences are that we now sum over all y G 2.^ 
and that arithmetic is modulo d. 

In the CSS formalism, a quantum code has a stabilizer 
group which is generated by a product of two distinct 
generating sets, the X-stabilizer generators, which are 
tensor products of X and I alone and the Z-stabilizer 
generators, which are tensor products of Z and I alone. 
In addition, one must define logical operators, and in the 
CSS formalism, logical encoded X operators consist of 
a tensor product of X and (optionally) I alone, and the 
logical encoded Z consists of a tensor product of Z and 
(optionally) I alone. Once the X-stabilizer generators 
and logical X operator (or operators if there are multiple 
encoded qubits) are defined, the remaining Z-stabilizer 
generators are fixed by conjugation relations which can be 
elegantly captured via the use of dual codes from classical 
coding theory. In the qudit setting, Gq will again define 
the X-stabilizer generators for the code, via the mapping 
k to X^, e.g. the row vector 000123 defines the operator 
70/0/0X(0 X^ ® X^, and Gi will similarly define 
the logical X operators for the code. 


B. Defining m*-orthogonality 

Bravyi and Haah m prove that a so-called triorthogo¬ 
nal code supports a transversal non-Clifford gate in the 
3rd level of the Clifford hierarchy. We now generalize 
Bravyi and Haah’s construction. First, we focus on gen¬ 
eralising from qubits to qudits. We defined qudit color 
codes using a star-bipartition, and it is useful to incor¬ 
porate that into the definition. We introduce the star 
sign-flip matrix F. Each column of matrix G corresponds 
to a qudit, and let us order these columns such that the 
first p columns correspond to unstarred qudits, and the 
latter n — p columns to starred qudits. We then define 
the n X n matrix F: 

E = diag( 1,..^.,1 ^.,- 1 ), (18) 

|t),| entries |t)*| entries 

where the first |t), | elements on the diagonal are 1 and 
the remaining elements on the diagonal are —1. We note 
that the ordering of the columns here is arbitrary, so we 
do not lose any generality by ordering the columns in 
this manner. The purpose of F will be to flip the sign of 
the entries of row vectors corresponding to the starred 
qudits, e.g. if 5 = {[ 5 ) 1 , [ 5 ) 2 ,. . ■, [g\n] is a row vector, we 
define g ■ F as the row vector where [gF]j = [g]j for the 
first p elements (corresponding to unstarred qudits) and 
[gF]j = —gj for the remaining elements (corresponding 
to starred vertices). We also define the weight | • • • | of a 
row vector as |g| = '^j[g]j. 

We now will define a significant generalisation of tri¬ 
orthogonal matrices, which we call TO*-orthogonality (pro¬ 
nounced “m star orthogonality”). 


(17) 
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Definition 2. An n x n matrix G over 1-d with a a 
bipartition of columns into is m*-orthogonal if 

both of the following conditions hold. 

1. The weight of every elementwise product of any m 
rows of GF (including repeated rows) is equal to 0, 
except for the following case: 

2. For each row of submatrix GiF, the weight of the 
row vector raised to the mth power elementwise is 

1 . 


where F is defined as in equation (181 and where the 
matrices’ columns are ordered to respect to the star- 
bipartition {z;,, ?;*}. 


Notice the above makes no mention of modular arith¬ 
metic. If we had instead defined the weights modulo 
d, we would have a weaker notion of m*-orthogonality. 
Color codes turn out to satisfy this stronger notion, and 
so this definition suffices for this paper. We note that 
the results below rely on the stronger form to deal with 
the exceptional cases, e.g. d = 3,6 for the T gate (see 
App.g. We also remark on the weaker notion to clar¬ 
ify that for d = 2 it exactly corresponds to the Brayvi- 
Haah definition of triorthogonality. Other applications 
of weak m*-orthogonality include quantum Reed-Muller 
codes [5115i]. 

A more symbolic statement of m*-orthogonality will 
prove useful in subsequent sections. We make use of of 
u o V symbol to denote element-wise products of vectors 
u and u, such that it has elements 


[9aO gb]j =[ga]j[9b]j, (19) 

which generalizes for an arbitrary number of vectors, e.g. 
for three vectors 


[9a o 9b o 9c] j = [9a] j [9b] j [9c] j ■ (20) 

Let us now consider a reformulation of points 1 and 2 of 
Def. For any set of m rows, we have a list of m row 
indices {a, 6,..., y}. m*-orthogonality demands that for 
all the 

Yll.9aF o gbF o ... o gyF]j = 0 (21) 

3 

unless all rows are identical and come from Gi, in which 
case the weight is unity. 


C. The m*-orthogonality of qudit color codes 

The concept of m*-orthogonality is a powerful tool for 
studying the properties of qudit color codes due to the 
following lemma 


To prove this, let us first prove a convenient lemma. 

Lemma 6. Any m*-orthogonal matrix G which includes, 
as one of its rows, the alTones vector is also m'*- 
orthogonal for all m' < m. 

Proof. This follows since an elementwise product of m 
vectors including the all-ones vector is equal to the ele¬ 
mentwise product of m — 1 vectors excluding the all-ones 
vector. □ 

Since the logical X for all color codes is the transversal 
X acting on all qudits, the matrix Gi for these codes 
contains the all-ones row. Hence to prove lemma we 
now only need to prove that a color code in m spatial 
dimensions is m*-orthogonal. 

Proof. Recall that the A-generators of the code are de¬ 
fined by yit'-cells. Taking as an example the 3D lattice 
in figure the A-stabilizer generators are defined by 
3-cells. Consider the following geometric properties of 
these cells. When q distinct /r'-cells intersect non-trivially, 
where they meet defines a cell of smaller dimension. For 
a general lattice, there is no further restriction on the 
dimension of this cell. However, for a /r-colex it is well 
known [T3i dsnni sg that the intersection of q objects 
of dimension fi' yields either an empty set or a cell of 
dimension p' — q -\- 1. For instance, the intersect of two 
neighboring /r'-cells defines a (/x' — l)-cell, where three 
such /x'-cells meet defines a (^' —2)-cell, and so on. Where 
fi' /x'-cells meet defines a 1-cell, or lattice edge. For exam¬ 
ple, in 3D, two adjacent 3-cells meet at a face, and three 
adjacent 3-cells meet at an edge (four adjacent 3-cells 
meet at a point). 

We use this geometric fact to prove the theorem. Any 
product of m row vectors in Go has a geometric repre¬ 
sentation as the intersection of the vertices of these cells. 
This corresponds to either an empty cell or cell of dimen¬ 
sion no less than p — m -\- 1. Each cell in the lattice, of 
any dimension greater than zero, has an equal number 
of starred and unstarred vertices (see Cor.[^. When we 
multiply the product row vector by F we invert the sign 
of the columns corresponding to the starred vertices, but 
this corresponds to precisely half of the non-zero elements 
in the vector. Hence the weight of the resultant vector 
is zero provided g, — m -\- 1 > 0. Therefore, the rows of 
Go satisfy the conditions for m*-orthogonality whenever 
m < pL. 

The other case to consider is the product where all 
m vectors in the product are the all-ones vector. The 
element-wise product results trivially in the all-ones vector. 
The number of unstarred qudits is one more than the 
starred qudits (see Cor.[^, hence the weight of this vector 
after multiplying by A is 1. Together with lemma this 
proves lemma □ 


Lemma 5. A qudit color code for any d > 2 defined 
on a lattice in spatial dimension fi where X-stabilizer 
generators are defined by g,'-cells with g' < g, and with a 
star-bipartition {u,,z;*} defined by the bipartition of the 
lattice, is a m*-orthogonal code for all m < g'. 


Note that the bounds in this lemma are tight. A code 
whose A-stabilizer generators are defined by /x'-cells is 
not {g' -\- l)*-orthogonal because (^+1) ^'-cells meet at a 
point, and a single vertex represents a vector which does 
not have zero weight. 
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V. TRANSVERSAL OPERATORS 


Our attention now turns to transversal gates. The qubit 
color codes are the only family of topological codes known 
to support transversal gates satisfying the Bravyi-Konig 
bound. We prove in this paper that (apart from a minority 
of special cases which need to be treated individually) the 
qudit codes we have defined also saturate this bound. 

Theorem 1. The qudit color codes on a fi-colex, where 
X-stabilizer generators are defined by g! < gi cells and 
where 0 (mod d), support transversal gates in the 

p'th level of the Clifford Hierarchy (CH) (which are not in 
the {p—l)th level of the CH), saturating the Bravyi-Konig 
bound. 

The qualification ^'! 7 ^ 0 (mod d) is due to lemma 
namely that when ff\ = 0 (mod d), gates of the form of 
equation are in the {ff — l)th level of the CH (and 
therefore does not saturate Bravyi-Konig) as discussed in 

Sec.lU 

For the most important p, = 3 case, however, we addi¬ 
tionally prove transversality for the exceptional d = 3 and 
d = 6 cases, and can thus state the unqualified result: 


Theorem 2. The qudit color codes on a p-colex, where 
X-stabilizer generators are defined by p' < p cells sup¬ 
port transversal non-Clifford gates in the 3rd level of the 
Clifford Hierarchy (CH) for all p' >3. 


We believe that the qualification in theorem can be re¬ 
moved, following a similar approach to proving theorem 
but we leave this for future work. 

To prove theorem we show that gates of the form of 
equation Q are transversal in all m*-orthogonal codes 
for polynomial degree r = m. To prove theorem ^ we 
also prove that gates of the form in equation ([7) are 
transversal in 3*-orthogonal codes. The theorems then 
follow by virtue of lemmas Bi and [Sl^gether with the 
known results for qubit color codes In particular, 

the transversal gates we consider are of star-conjugate 
transversal form, see equation El 


Lemma 7. A m*-orthogonal code has a transversal im¬ 
plementation of the unitary gate Rf^ such that: 


Rfrr,= ( 22 ) 

implements Rf^ on the code space. 

We prove this lemma in App. As a simplified presen¬ 
tation for the main text of the paper, we prove the special 
and important case of r = 3. The general proof follows 
the same approach, but is a little notationally unwieldy. 

Proof. 

Lemma 8. A 3*-orthogonal code has a transversal im¬ 
plementation of the unitary gate T = Rj^ such that: 


f = T[v.](^T*[v4 (23) 


implements T on the code space. 

We now examine the conditions that the star-conjugate 
transversal T implements a logical T. For this to be 
true, the phases for each computational basis component 
of each logical codeword in equation 0 ’ must agree 
with the phases defining the gate in equation §. Noting 
that the same phase applies to all terms in the sum in 
equation ( [l7| ), we can write 

f\{x,yf ■G)=u;-"\{x,yf -G), (24) 


where x € J-d and y € where s is the number of 
stabilizer generators. In other words we fix the number 
of logical qubits k = 1. 

Applying equations Q and (23) to the state on the 
right hand side of equation (241, we recover the following 
relationship between the phases on both sides of the 
equation: 


((x, yf' ■G-F)o {{x, y)'^ -G- F)o 
((x, yY' ■ G ■ F) = x^ (mod d). 


(25) 


where we have right multiplied the F matrix to vector 
{x, yY ■ G. It is useful to define z = {x, y) so 


Y,{{z'^ ■G-F)o{z^ ■G-F)o{z^-G-FY 

3 

= x^ (mod d). (26) 

Now define z = ZoCa where Zq is the value of the 
ath position of the vector z, and where Ca are the basis 
vectors over ■ For example ei = (1,0,... ,0), 62 = 
(0,1,..., 0) and so on. Re-writing the left hand side of 
equation ( [2^ we obtain 

■ G ■ F)o ■ G ■ F)o ■ G ■ F)), 

■ G ■ ^ (mod d). 

Note that • C? is the ath row of the matrix G, which 
we denote ga. Therefore 


^((z^.G-K)o(z^.G-F)o(z^.G-F)), 




(mod d). 
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we can re-write this as 

3 

s+1 

= X! X! ( (^“5“ ■ ° ■ F) O 

j a,b,c—l 

{zc9c-F))^ (mod d), 

s-l-l 

= E E [{gaO gbO 9c) ■ F]-ZaZbZc (modd). (27) 

j a^b,c—l 

Summing over j gives 

■G-F)o{z^ ■G-F)o{z^-G-F)), 

3 

s+1 

= E \{9aO gbO gc) ■ F\zaZbZc (modd). (28) 

a,b,c=l 

We break the sum into two parts, so that 

^((^^ • G • F) o (z^ . G • 1^) o (z^ . G • F)), 

3 

= E KdaOgaOga) ■ F\zl + 

a 

\{gaO gbO gc) ■ F\ZaZbZc 

{a,b,c\—i{a—b—c)} 

(mod d). (29) 

Now we use the definition of m*-orthogonal matrices. 
This ensures that the first term is equal to (mod d) 
and the second term is equal to zero. Hence equation ( |2^ 
is satisfied, which completes the proof of lemma □ 

For larger m, the proof follows the same lines but is with 
more cumbersome algebra. The proof of lemma which 
we present in App. A] completes the proof of theorem 
To prove theorem [2 we require one further lemma: 

Lemma 9. A 3*-orthogonal qudit code where d = 3 or 
d = 6 has a transversal implementation of the unitary 
gate T 3 g defined in equation 0 such that: 

T3,6 = T3,6[z;.](8)T3*gK] (30) 

The proof is presented in App. completing the proof 
of theorem]^ Finally, we remark, that as CSS codes, all 
color codes admit a transversal A(A). We conclude this 
section with a table of star-conjugate transversal gates 
in qudit color codes of different spatial dimension, see 
table m 

VI. GAUGE FIXING TO IMPLEMENT THE 
LOGIGAL HADAMARD 

In this section and Sec. EIIl we show how a universal 
gate set can be achieved via gauge fixing, or by defin¬ 
ing a gauge color code. Although these results follow 


directly from the qubit cases presented in [ 2 S] and m 
without the need for any further technical innovations, we 
include a discussion and explanation of these techniques 
for completeness. 

First we consider gauge fixing. We shall define a sub¬ 
system code occupying the same lattice as the color code. 
This code is capable of realising the logical Hadamard 
gate on the same encoded qudit as the color code, while a 
set of gauge qudits become corrupted. Fortunately there 
is a simple protocol to fix the corrupted gauge qudits to 
re-initialize the color code. 


In a [[n, fc,d]] subsystem code [55H5S] there are three 
sets of operators that define the code. First, a set of gauge 
generators is defined, a set of Pauli operators defined on 
the lattice that do not necessarily mutually commute. 
These generate the gauge group A. The center of the 
gauge group (the elements of the group which commute 
with all other elements) represents the stabilizer of the 
code S (up to a free choice of plus or minus phases which 
must be chosen to define the zero-error code space). The 
final set of operators to define is the set of logical Pauli 
operators on the code space, the Pauli operators that 
commute with A. Subsystem codes can be understood as 
stabilizer codes in which some of the logical qudits have 
been demoted to gauge qubits. These gauge qubits are 
not used to carry information and can be corrupted or 
measured during operations on the code. 

The logical operators X and Z are products of X and 
Z (and X* and Z*) operators acting on every qudit. 
Therefore, a natural candidate for transversal 77 is a 
global star-conjugate Hadamard. However, for color codes 
of // > 2 it can be seen from the structure of the stabilizer 
generators that such a unitary will not leave the code 
space invariant. In particular, the X stabilizers act on 
the /x'-cells of the lattice whereas the Z stabilizers act 
on the {pL — ji' -\- 2)-cells of the lattice. Therefore unless 
g' = pL — g' + 2, the Z stabilizers cannot be transformed 
into valid X (cell) stabilizers, by a global Hadamard. The 
solution to this is to switch between the stabilizer color 
codes and a different code in which X stabilizer generators 
and Z stabilizer generators have identical support. 


In g spatial dimensions the subsystem code stabilizers 
are defined as the star-conjugate X and Z operators acting 
on the qudits contained in a /x-cell of the lattice. The 
logical operators matching the color code logical operators 
are the transversal X and Z defined in equations (141. In 
addition there are X and Z gauge operators defined on 
the {g — l)-cells, {g — 2 )-cells, and so on to the 2 -cells of 
the lattice, such that they commute with the stabilizer 
group. They can be identified in w-commuting pairs, with 
each pair protecting a gauge qudit. 


In higher spatial dimensions there are a larger number 
of gauge generators and hence more gauge qudits than in 
lower spatial dimensions. Nevertheless, the basic protocol 
for realising the logical Hadamard gate in the subsystem 
code and fixing the gauge to return to the stabilizer code 
is independent of the spatial dimensions. 
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Logical operator 

Defined in Eqnation 

Star-conjngate transversal implementation 

Spatial dimension 

H 

A(A) 

S 

T 

Rfr 

eqnation ( 
eqnation ( 
eqnation ( 
eqnation ( 
eqnation ( 

3 

5 

4 

6 
8 

1 

1 

1 

1 

1 

H[v,] (g) 

Applied blockwise transversally 

S[n.] (g) S'*]!!*] 

T[v,] (g) r*[ii*] 

Rffiv,] (g) 

2 

> 2 

> 2 

> 3 

> r 


TABLE I: A table illustrating the star-conjugate transversal implementation of some important logical gates and the spatial 
dimensions of the codes which support them. 


The star-conjugated transversal Hadamard gate, 

H = H = H[v,]®H*[v.], (31) 

transforms the transversal logical operators and the stabi¬ 
lizers of the code correctly. However the gauge generators 
are corrupted under the action of H, meaning that Zj 
(the Z logical operator for qudit j) is transformed not to 
Xj as desired but instead to the operator for a different 
qudit, Xk, causing the gauge qudits to become entangled. 

These corrupted gauge operators must be fixed, to 
return the code to the original stabilizer code in which 
the other transversal operators can be realized. This hxing 
is enacted by measuring the stabilizer generators. The 
measured eigenvalues will indicate the Pauli correction 
operator that should be applied—this is an operator that 
restores the corrupted gauge qudit to its initial state while 
commuting with the other logical operators. 

Once the corrections are applied, the gauge qudits are 
re-initialized to their original state so that any entan¬ 
glement between the gauge qudits is destroyed. This 
is equivalent to the stabilizer color code with a logical 
Hadamard applied to the encoded qudit. 


VII. QUDIT GAUGE COLOR CODES 

Another approach to the implementation of the 
transversal Hadamard is to treat the subsystem code 
as the base code, and to achieve the transversal non- 
Clifford gates by gauge fixing to the (stabilizer) color 
code. This switch in perspective leads to the gauge color 
codes, recently introduced by Bombin [15] . Far more than 
a change in perspective, however, gauge color codes gain 
many useful new properties, which the original color codes 
do not possess. 

The first advantage of the gauge color code construc¬ 
tion is that the outcome of stabilizer measurements can 
be inferred from measurements of the gauge operators. 
This offers a practical advantage, since the weight of the 
gauge generators can be significantly smaller than sta¬ 
bilizer generators. For example, in 3D the X stabilizer 
generators would include 24-body operators associated 
with truncated octahedrons, discussed in Sec. Imcl re¬ 
quiring coherent 24-body measurement. In the equivalent 
gauge color code, however, the X stabilizers can be decom¬ 
posed into 4- or 6-body X gauge operator measurements, 
a much more practically feasible proposition. 


A second advantage (in the ^ = 3 case) is that the 
outcome of stabilizer measurements is more robustly en¬ 
coded in the gauge generator measurement. The extra 
redundancy leads to to “single-shot fault-tolerant error 
correction” m — meaning that measurement errors can 
be accounted for without the need for repeated measure¬ 
ments. 

A third advantage is that one can use gauge fixing 
to fault tolerantly convert m between codes of differ¬ 
ent spatial dimension [36], allowing one to combine the 
advantageous features of both. For example, one could 
perform all Clifford computations in a 2-D code while 
only resorting to a 3-D code to achieve a non-Clifford 
gate. 

We can define qudit gauge color codes in any spatial 
dimension analogously to Bombin’s qubit codes |15| . The 
only difference is that operators are star-conjugated ac¬ 
cording to the star-bipartition of the lattice. 

The most studied gauge color code exists on a 3D lattice 
and we shall describe the qudit analogue as an example. 
Gauge generators consist of star-conjugated face operators 
of both X and Z-type. The stabilizer is generated by 3-cell 
operators whose values can be obtained from the products 
of the face operators. The star-conjugate transversality of 
the Hadamard gate for this code is trivial via its definition. 
To achieve a transversal T gate one must gauge fix back 
to the stabilizer color code by promoting Z face operators 
to the stabilizer. 

Although the basic properties of color codes transfer 
immediately to the general qudit case via this dehnition 
and by the results presented in previous sections in this 
paper, the performance of single-shot error correction is 
less straightforward to analyze, and we leave a detailed 
study to further work. Single-shot error correction works, 
in the qubit case, because the relationship between gauge 
and stabilizer measurement outcomes in the code is such 
that the values of the stabilizer are mapped to delocal¬ 
ized structures (branching points in string-like extended 
structures) in the gauge outcomes. This arises because 
the parity of the outcomes of the face operators around 
the cell must be equal to the eigenvalue of the stabilizer 
representing the cell, even parity in for a trivial syndrome, 
and odd-parity for the non-trivial case. 

In qudit codes, the stabilizer measurement outcome 
may take any of d — 1 non-trivial values. The outcomes 
of gauge measurements must sum to these values. Thus 
the relationship between gauge measurements and stabi¬ 
lizer measurements is more complicated than the simple 
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FIG. 4: (Color online) A cell of the 3D gauge color code. 
Its faces are 3-colorable and colored red, green and blue in 
this example. The stabilizer outcomes for the cell can be 
constructed by adding together (mod d) the outcomes of the 
gauge measurements for each of the red faces, or each of the 
blue faces, or each of the green faces. 

“branching point or no branching point” behavior seen in 
the qubit code. This extra structure, however, should 
represent additional information about the location and 
charge of stabilizer outcomes which a decoder could ex¬ 
ploit. Thus, with the development of a suitable decoder, 
we expect single-shot error correction robust to measure¬ 
ment error to be achieved in qudit gauge color codes, 
but we leave its the construction of such a decoder, and 
analysis of its performance, to future work. 

VIII. DISCUSSION 

We have shown how any existing qubit color code can 
be generalized to support a qudit color code in arbitrary 
spatial dimensions, and with an arbitrarily large distance. 
To do this we introduced the notion of a star-conjugate 
lattice and showed how this construction allows the Clif¬ 
ford phase, controlled-V and non-Clifford phase gates to 
be implemented transversally. The set of Clifford group 
generators is completed by the Hadamard, which can 
be implemented in a subsystem code. Techniques for 
switching between the two codes have been outlined. An 
additional advantage of considering the color codes as 
gauge codes is reduction of the weight of measurements 
required to obtain the error syndrome. 

There are still many open questions remaining—here 
we shall outline some of the ones we find most interesting. 
Although a fault-tolerant decoder has been implemented 


for a 2D code, no decoder has, to our knowledge, been 
proposed for color codes in higher spatial dimensions. In 
Ref. [37] a decoder for a 3D gauge color code is discussed, 
unfortunately however, a classical algorithm capable of 
performing the decoding is not specified. One would 
expect that in a fully fault-tolerant simulation the gate 
error threshold for the ^-dimensional color codes would 
be prohibitively low, since the stabilizers defined on the 
/i-cells would require circuits of many time steps to mea¬ 
sure. There is hope however that by switching to the 
corresponding gauge color code the measurements require 
fewer gates and therefore the threshold may recover m- 
The development of efficient decoders for qudit gauge 
color codes is therefore of significant interest. 

The construction of the transversal gates in the color 
codes and gauge color codes presented in this work holds 
for qudits of any dimension, d. However, for non-prime d it 
is not clear that S, H and A(A) are sufficient to generate 
the Clifford group. We desire a better understanding of 
universality and the Clifford hierarchy when using systems 
of non-prime dimension. 

We have identified some other lattices that can sup¬ 
port a color code in 3D. It would be interesting to see 
a comprehensive list of possible lattices. Furthermore, 
our star-conjugate construction is valid for all of these 
lattices, but there may be one that offers advantages in 
terms of the necessary qudit overhead to achieve a desired 
distance. 

Finally, there are many Reed-Muller codes known for 
qudits for which no equivalent color code is known jS] [33] . 
Is it possible to use one of these to construct a topological 
code similar to a color code and study its properties? Here, 
we found codes where the distance is purely topological 
in nature, and so depends on the lattice not on whether 
we use qudits. In contrast, qudit Reed-Muller codes 
provide an improved distance that is algebraic rather 
than topological in origin. The potential exists for qudit 
quantum codes which abstract topological ideas to a more 
general setting to generate novel codes of arbitrary code 
distance with a rich family of transversal logical gates. 
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Appendix A: Proof of lemma 


Lemma 7. A m*-orthogonal code has a transversal im¬ 
plementation of the unitary gate Rf^ such that: 


= Rfrr. [f.] Rf^ [w] 
implements Rf^ on the code space. 


(Al) 


Proof. For equation (Al) to be true, for all z = {x, y) we 
have 


■ G) = ■ G) ■ 


(A2) 


Applying equations ^ and (22) to the state on the right 
hand side of equation (A2), the phases on both sides of 
the equation are equal if 

(z^ • G • F)°^ = (mod d), (A3) 

where the notation (.)°™ means a circle-product (see equa¬ 
tion (iHl)) taken between m identical elements. Proceeding 
using the same notation as in Sec. |V]we find 


^((z^.G-F)°™), 

j 


(mod d), 


= II| H ^sga.F 


(mod d), (A4) 


where [a,..., d]"” is a vector of length c with elements that 
take values between a and b, and where gs is the bitwise 
product of m rows of G, so 


gs = gai o ga2 ° ° ga 


(A5) 
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We can re-write this as 

-G-FT'-), 

3 

= X! X! [ 9 a-F\jZs (modd). (A 6 ) 

3 aG [l, ■ ■ ■ 

Summing over j gives 

Y,{(.z^-G-Fr^), 

3 

= X! \9s-F\zs (modd). (A7) 

a£ [1 ,... ,/l+s]”^ 


where 7 ^ = a; and where the function in the exponent 
is evaluated in regular arithmetic (or equivalently modulo 
3d). We shall prove lemma for d = 3 and d = 6 
separately, since each requires calculation in arithmetic 
with a different modularity. 


1. d = 3 

Proof. First we consider the case of qutrits, where d = 3. 
In this case, 7 = e*“5“. Hence, our calculation will make 
use of addition modulo 9. Proceeding as before we start 
with: 


We now use the definition of a m*-orthogonal code. 
This means that all terms satisfying ^(oi = ... = am) 
in the summation on the right hand side vanish. This 
implies that the only non-zero contribution is represented 
by d = (a,a,...,a), a m-element vector where every 
element is a. Recall also that such terms vanish unless 
a = 1 , in other words it is the m-fold elementwise product 
of the vector in Gi. We note that zg = zff, leaving only 


Y,{{z^-G-F)°n, = \9s-F\zT (modd), 

3 

= (mod d), (A 8 ) 


which proves the desired equation (A3) is satisfied. □ 


Appendix B: Proof of lemma 

Lemma 9. A 3*-orthogonal qudit code where d = 3 or 
d = 6 has a transversal implementation of the unitary 
gate defined in equation ([^ such that: 


fa,6 • G) = 7 "' • G). (B3) 

The phases on both sides of this equation match if 

^((z^ • G • F)°^)j = x'^ (mod 9). (B4) 

3 

This expression is a mix of modulo 3 and modulo 9 arith¬ 
metic so we must proceed with care. The matrix multi¬ 
plication is mod 3, whereas all other arithmetic is mod 9. 
Re-writing the left hand side as 

■ G ■ F)°% (B5) 

3 

s+1 

= H {[9a-F]jZa (mod3)}o 

j a,b,c—l 

{[gb-F]jZb (mod3)}o 

{[gc ■ F]jZc (mod 3)} , (B 6 ) 

we can use the following identity to express the modulo 3 
reduction in terms of standard addition 


^ 3,6 = 73.6b»] ® T'a.ef* 


where 


f3.6 = 1-?'^ 

je'^d 


(Bl) 


a (mod n) = a — I — In. 

n 


(B7) 


Setting aside the modulo 9 arithmetic for the moment 
and evaluating this expression first in standard arithmetic, 
(B2) the left hand side now expands to the following form: 
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s+l 


J2iiz^-G-Fr)j = j2 E 


j a,fc,c=l 

o (idb ■ F]jZb - 3[ 
s+l 


o L, ■ F],z, - si ^9c-n^c 


= E E [(5a O 5& O 9c) ■ F]jZaZbZc - 9[(ga O gb) ' F]jZaZbl^^^^^^^\ 


j a,b,c—l 


+ 27[ga ■ 

_ 2y^ ’ ^]i^a J ^\3bj_F]j_^^ ^ [9c • F]jZc ^ 


r 


The last three terms in this summation are integer We expand the right hand side and find that 
multiples of 9 and are thus equal 0 (mod 9). Thus if we 
reimpose modulo 9 arithmetic we are left with: 


■ G ■ F)°^), (mod 9) 


s+l 

= E E °9b° 9c) ■ F]jZaZbZc 

j a^h,c—l 

(mod 9) (B 8 ) 

Finally, invoking TO*-orthogonality completes the proof. 

□ 

Note that we are using the fact that the definition 
of TO*-orthogonalityis stated in regular arithmetic which 
implies the modulo 9 orthogonality used here. 


E((2^ • G • F)°% 


s+l 


= E E [(5a o 5& o 9c) ■ F]jZaZbZc (mod 18). 

j a,b,c—l 

(BIO) 


2 . d = > 


Proof. When d = 6 the phase in is 7 = e* is and thus 
proceed with modulo 18 arithmetic. 


T3,6 \z^ • G) = 7 " \z^ • G) . 


Converting modulo 18 arithmetic to normal arithmetic 
(B9) we find, 


J2iiz^.G-F)°%=Y: E [(5a O 56 O 9c) ■ F]jZaZbZc - 18[(5a O db) ' F]j ZaZbl^^^^-^li^ ^ 

j j a,b,c=l 

+ 108[5a • l ^3c-Fhz, ^ 

— 216[J^(hjEk^j ' ^]j^6 j 1 ^ [5c • F]jZc ^ 


Casting this back into modulo 18 arithmetic we can cancel the terms equal to integer multiples of 18, leaving: 

^(( 0 ^ • G • F)°% (mod 18) 


s+l 


= E E [(5a o 56 o 9c) ■ F]jZaZbZc (mod 18) 

j a,b,C—l 


(BID 
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Finally, via m*-ortliogonality we recover 

■ G ■ (modl8)=a;3 (mod 18) (B12) 

3 

completing the proof. □ 



